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Abstract 

We compare within an unifying formalism the dynamical properties of modulated 
and composite aperiodic (incommensurate) crystals. We discuss the concept of inner 
polarization and we define an inner polarization parameter (5 that distinguishes between 
different acoustic modes of aperiodic crystals. Although this concept has its limitations, 
we show that it can be used to extract valuable information from neutron coherent 
inelastic scattering experiments. Within certain conditions, the ratio between the 
dynamic and the static structure factors at various Bragg peaks depends only on f3. 
We show how the knowledge of (3 for modes of an unknown structure can be used to 
decide whether the structure is composite or modulated. The same information can be 
used to predict scattered intensity within unexplored regions of the reciprocal space, 
being thus a guide for experiments. 



1 Introduction 

Aperiodic crystals are long-range ordered structures whose diffraction patterns are made of 
Bragg peaks. The difference with respect to normal crystals is that they are not periodic 
and one needs n = 3 + D (D is the dimension of the internal space) basis vectors to index 
the positions of peaks in reciprocal space: 
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q = y^»q», Zi e z (l) 

8=1 

This property has been exploited in crystallography by embedding the structure in a 
superspace crystal of dimension n. Atomic positions are the intersections between the D- 
dimensional atomic surfaces and the 3-dimensional physical space. For periodic crystals, 
the Floquet-Bloch theory reduces the 3N dimensional eigenvector problem of lattice modes 
(N is the number of atoms) to a 3N U dimensional problem (N u being the number of atoms 
inside the unit cell). The similar approach can be used for aperiodic crystals, but in this 
case the number of equations is not reduced; because of incommensurability N u = N. Nev- 
ertheless, as we show in this paper, superspace concepts lead to important simplifications 
in the description of modes, that can be used in the analysis of inelastic neutron scattering 
experiments. 

At present, three classes of aperiodic crystals are known: incommensurate modulated 
crystals, incommensurate composites and quasicrystals @> H> S H- Atomic surfaces are 
discontinuous for quasicrystals, but they can be continuous for modulated crystals and for 
composites. The borders between the three classes are not so clear as one may think. The 
classification difficulties between modulated and composites crystals are notorious. If dis- 
continuous atomic surfaces and occupational modulations are allowed ^ any modulated 
structure can be seen as a composite one and vice-versa. The choice can be particularly 
difficult in compounds with a complex structure. For instance, recent structural neutron 
investigations of Bi-2212 superconductors J7| emphasize the difficulty of deciding whether 
these materials should be called composite or modulated. A natural question arises. Can 
dynamics as probed by inelastic neutron scattering experiments [ |TU| , [| give an answer to 
this classification problem? 

In order to answer this question we need to know the typical dynamical responses of 
different aperiodic crystals. Several simplified model systems will be investigated, hoping 
that more complex systems behave in resemblance to one of these. Of course, it is also 
possible that compounds that structurally are at the border between the three classes of 
aperiodic crystals have also special dynamical properties, at the border between the ones 
studied here, in which case the classification problem could be really intractable. 

In this paper, we discuss the linear (small amplitude) excitations of displacively modu- 
lated and composite crystals. The case of occupational modulations and that of non-linear 
(large amplitude) excitations will be studied elsewhere. 

We focus our discussion on the dynamics of low frequency excitations. The n-dimensional 
superspace picture suggests that aperiodic crystals have n branches of propagating hydro- 
dynamic modes. Only three of those are true Goldstone modes, associated to the broken 
translation symmetry of the Hamiltonian. These three branches are the general acoustic 
phonons, present in any solid. The other p branches are called phasons because they are 
associated to the broken "phase" symmetry of the ground states manifold (which is not a 
symmetry of the Hamiltonian). 

The hydrodynamics of composite aperiodic crystals has already been the object of several 
theoretical papers [L2|, [13], Although they offer rather complete results (including 
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damping) these papers suffer from their generality and from the use of various phenomenolog- 
ical parameters. It is difficult to extract from them the consequences of incommensurability 
on dynamics and to find the answers to simple (yet important) questions such as which is 
the contribution of each mode to the neutron inelastic scattered intensity and which are the 
regions in the reciprocal space where these contributions are most intense. Modes observ- 
ability is not only determined by the (integrated) intensity of the peaks, but also by the line 
shape. The hotly debated problem of damping, whether intrinsic as mode-mode coupling 
[ |3^ or extrinsic via coupling to defects, is left out from our discussion. 

The superspace picture suggests that there is a tight connection between structure and 
dynamics. From the very beginning of our paper we exploit this connection and our purpose 
is to find a small set of quantities that can be easily interpreted theoretically and measured 
experimentally, and that facilitate the classification of aperiodic crystals from the point of 
view of their dynamical properties. 

2 Ground states, modulation functions and low fre- 
quency excitations 

Throughout this paper we shall deal with ID models, otherwise we shall stay as general as 
possible. 

The positions of the atoms, corresponding to a ground state of a displacively modulated 
incommensurate crystal, read ||: 



y n , g = nai + 5 g + f g {nai + 8 g ) 
where f g (y + a 2 ) = f g (y) 

The indexes n G Z, g — 1, n\) are for the unit cells of the basic (non-modulated) structure 
and for the atoms within this unit cell, respectively. f g are periodic modulation functions of 
period a 2 , &\ is the period of the basic structure, the ratio a = a 2 /a± is irrational. 

A composite has several subsystems. We consider here the simplest case when there are 
two subsystems. The positions of atoms belonging to the two subsystems read [I5|: 



y ( Z = ma 2 + ^ ) + fl 2 \ma 2 + 5^) 



(3) 

and fh'(y + a = fh'(y)> h = i,n 2 



where fW(y + a 2 ) = ffl(y), g=\,n x 



"(2)/ , . „ X _ ,(2), 



The modulation function of one subsystem has the periodicity that the other subsystem 
had before modulation. 

The ground states of incommensurate modulated and composite crystals are degenerated. 
This can be explained by the existence of symmetry groups that leave the configuration 
energy invariant fl5| : 
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• Physical space translations, that are uniform displacements. 

T\(y n ,g) = Vn,g + A (4) 

where A G R 

• Discrete inner space translations, that are combinations of discrete uniform shifts and 
relabelling of atomic positions. 

Pr,s{yn,g) = Un+r,g ~ Ta X (6) 

PrAVnii Vm,h) = (yn+r,g ~ ™1, ~ ™l) ( 7 ) 

where [] r, s G Z. 

The transformations P T)S change the phase of the modulation functions and (only for 
composites) produce a relative uniform displacement of the subsystems: 

P r , s ({nai + S g + f g {nai + S g )}) = 
{n ai + 5 g + f g {na x + 5 g + 5')} (8) 

P r ,,({na 1 + 4 1) + /, (1) (na 1 + 4 1) ), 
na 2 + 5^ + ft\na 2 + 5^)}) = 
{na 1 + 4 1) + tf ) (na 1 + 4 1) +(5'), 
na 2 + 5f - 5' + /f (na 2 + < 2 } - 5') } (9) 

where the phase variation is: 



5' = ra x — sa 2 (10) 
Physical space and inner space translations have a simple geometrical interpretation in 



a superspace embedding [0 where they represent uniform translations of the superspace 
crystal. Physical space translations are parallel to the physical space. Inner space transla- 
tions make a non-zero angle 9 with the physical space (Fig.|l|). The physical effect of phase 
translations depends on 9 but also on 71,72, which are the angles made by basic superspace 
lattice directions and the physical space. In Fig.|I] these lattice directions define the fixed 
repeat distances a\ and a 2 along the physical space, while 71,72 are free to change. It is 
convenient to chose 71 = tt/2. In this case inner space translations defined by Eqs. |^|] 
correspond to 9 = tt/2 (they are orthogonal to the physical space). 

1 the action of s is trivial for modulated crystals 
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More generally, we may define composed translations as combinations TxP r , s that cor- 
respond to translations in superspace along directions making an angle 9 with the physical 
space. If 7x = tt/2 then one has: 



tan(9) _ 5' 

tan(<y 2 ) A 1 ' 

Because 02/01 is irrational, the phase variations 5' (EqJT0|) can approximate with arbitrary 
precision any real number (they form a dense set). If the modulation functions are smooth 
(the so-called analytic regime |18j) the group of inner space translations can be extended 
to a continuous group {-FWeR just by using Eqs.§|5] previously written for 5 f = ra x — 
sa 2 , r, s G Z and extend them by continuity to all real values of 5'. Physical space and inner 
space translations transform ground states into ground states. As a consequence of the 2D 
continuous degeneracy of the ground states manifold there are two hydrodynamic modes. In 
the limit of infinite wavelength k — > 0, the displacements involved by these modes are the 
infinitesimal physical space and inner space translations 



Infinitesimal physical space translations 

u n,g = e (12) 

,,gi u rn,hi 

Infinitesimal inner space translations 



{«£U£j = foe} (13) 



u n , g = V ^(na 1 + 6 g ) (14) 
iJX) J?) \ - 

\ a n,gi U> mh f — 

{V^ina, + #>), - V [l + f(ma 2 + 8® )]} (15) 

Let us define the inner polarization as a vector in superspace: P = (e, 77 tan(7 2 )) . A 
vibration of the superspace crystal in the direction P will correspond to a combination of 
infinitesimal physical space and inner space translations with coefficients e,r] (Figjl]). The 
inner polarization should not be understood as the direction of atom displacements (this is 
the physical polarization, always contained in the physical space), but as a vector belonging to 
the superspace (abstract construction describing internal degrees of freedom in a geometrical 
way). When the angle 9 between the inner polarization and the physical space changes, it 
is only the sequence of signs and magnitudes of atom displacements that changes, their 
direction remaining the same (Fig.0). 

The relation between 9, e, r\ is analogous to Eq.|Tl"|: 

tan{9) = V _ 
tan (72) e 

Let us define the polarization parameter (5 as: 
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= ^ (17) 

2e — 77 2tan(7 2 ) -tan(0) v ; 

In composites (3 may be called "sliding parameter" because at k = 0,u = this is 
proportional to the relative displacement between the centers of mass of the subsystems. 

p = —£ w— ( 18 ) 

< ugj > + < uZ\ > 

the averages < * > being with respect to all atoms in a subsystem. 

From Eq.[l8] it follows that in composites the inner polarization parameter and the par- 
ticipation ratios TZ\,TZ2 (that represent the fraction of the total "energy" concentrated on 
each subsystem) are related. In the zero th order of the modulation amplitude one has: 



a < \u„'g\ > a 



,(1)|2 

K, = - - ' ' gl - « - (19) 

«<kSi 2 > + <kSj 2 > a+t2 



where t = ttz, ol = a 2 /ai. 



1+/3' 

The modes are concentrated on the first subsystem if TZi = 1, IZ2 = 1 — TZi = 0. This 
occurs when /3 — 1, 9 — 72, thus when the inner polarization is along the atomic surfaces of 
the second subsystem. When (3 = — 1, 6 = n/2 (inner polarization along the atomic surfaces 
of the first subsystem) IZi = 0, 7?.2 = 1, hence the modes are concentrated on the second 
subsystem. 

In the above reasoning we supposed that modulation functions are smooth and that 
u = 0. We would like to know if P and (3 can be defined also elsewhere. 

Using group theory reasonings one may show that modes in aperiodic crystals are 
generalized Bloch waves of the type: 



V"n,g,K,y 

= exp[iK(nai + dg^Ug^in^ - y/tan^)} (20) 
for modulated crystals 
\u {1) u {2) \ - 

\ a n,g,K,yi u m,h,K,yS ~ 

%[nai — y/ tan(7 2 )], 

exp[iK(ma 2 + 5^)]U^ K [ma 2 + y/taa( T2 )]} (21) 
for composites 

K is the reduced wave vector (K = k — 2-nir ja\ + 5/02)), y is the internal space coordinate 
(orthogonal to the physical space), 72 is the angle in Fig.|l] (and 71 = vr/2). For each value 
of y one has an equivalent realization of the aperiodic crystal, so one may take y = in 
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order to obtain the actual displacements. U 9i k, uj^ K are periodic functions of y, having 
periods 02,02 and a±, respectively. We shall call them hull functions and it is useful to 
picture them in superspace. At K = 0, the hull functions play in dynamics the same role as 
the modulation functions play in statics, i.e. their intersection with the physical space give 
the displacements of atoms (Fig|J). At K 7^ the hull function displacements are multiplied 
by sinusoidal plane waves of wavelength 2tt/K (Eqs. p0"|j21~|) . Different possible hull functions 
are represented for modulated crystals in Fig.Q. In composites, one should imagine two 
sets of atomic surfaces and equal polarization vectors P for the two. As can be understood 
from FigfJ, the inner polarization can be defined if and only if the atomic surfaces are not 
distorted. The atomic surfaces are undistorted in a vibration mode if the motion can be 
described as a rigid motion in superspace. In physical space this comes down to a motion 
that does not change the local isomorphism class (the set of local atomic configurations 
remains the same). In the analytic regimes this implies that the hull functions are of the 
following form: 



Ug!K (y) = 6(K)+ V (K)^(y + 6 g ) (22) 



U^Uy) = e(K) + v(K) d ^[y + #>] (23) 



Uj?l(y) = e(K) - rj(K) - v(K) d ^[y + of } ] (24) 



Eqs. p2||2~3"| , p4| are exact for K = 0, to = in the analytic regime. We conjecture that they 



are fulfilled with good accuracy for small K, co. 

In Sec.|] we shall check the validity of the undistorted atomic surfaces hypothesis for 
different models of composite and modulated crystals. For the analysis of these models we 
shall also need a set of practical formulas to compute parameters T](K),e(K), (3(K). If Eqs. 
2^j2^j2^ are valid then one has: 



e{K) =< U g>K {nax) > (25) 

= < u a,K( na ^p( na i+ s 9> ^6) 



(3{K) 



<(S) 2 > 



2e(K)- V (K) 
<C/ 9 ,K(n«i)^(nai+<5 9 )> 



2 <(S-) 2 ><' 7 9.^( na i)>-<^ 9 .A'(nai)|i-(nai+ ( 5 9 )> 

for modulated crystals 
e(K) =< U { g % > (28) 

V(K) =< U<% >-< Ui% > (29) 

R(K\ — _ <j7 i-ff > ~ < ^,^ > (Qf\\ 

for composite crystals 
Several things are worth to be noticed: 
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After replacing the displacements by the hull function displacements, i.e. after elimi- 
nating the sinusoidal plane wave of wavelength 2tt/K, e(K) , r)(K) , (3(K) are calculated 
in the same way as they were at K = (compare Eqs . p~6| , p~T| , p~8| and Eqs. [27|j30|) . 

In modulated crystals, e(K) is the displacement of the centre of mass and r](k) is the 
linear regression coefficient between the hull function displacements and the derivative 
of the modulation function. 

In composites, e(K) is the displacement of the centre of mass of one subsystem and 
r](K) is the relative displacement between the centres of mass of the two subsystems. 

Contrary to the physical polarization which does not change inside a branch, the inner 
polarization expresses phonon/phason coupling and may depend on K. 

Generally, Eqs. |22],|23],|24| are only approximate (atomic surfaces are distorted) and their 
sum of squares errors [EH] are given by: 



SSE(K) = J2 n jU g>K (n ai ) - e(K)- 

- v (K) d £(na 1 + 5 g )\* (31) 
for modulated crystals 
SSEW = ^\U®{na 1 )-e(K)- 

- V (K) d 4^(na 1 + 5 9 1) )\ 2 + 

V (K) d ^-(ma 2 + 5^ (32) 
for composite crystals 

Eqs.|25|, |26|, |28|, p9| provide least squares regression coefficients corresponding to mini- 



mum SSE (Eqs.jJl]^). The inner polarization parameter follows from Eqs.p7|,30 



Traditionally, modes with 6 = are called acoustical phonons. The name emphasizes 
the fact that these modes are generic Goldstone modes, occurring in any solid. Modes with 
6 7^ are usually called phasons, because they involve non-zero phase fluctuations of the 
modulations. We prefer to call phason only the mode whose hull function conserves the 
center of mass of the solid. In modulated crystals the phason has e = 0,(3 = —1,6 = it/2 
corresponding to a vibration of the superspace crystal orthogonal to the physical space, hence 
along the mean atomic surfaces. The atomic displacements have zero average and conserve 
the centre of mass. In composites the phason has (3 = (p^ + p^')J{p^ 2 ' — p^),tan{6) = 
(1 + j ' p^)tan{^2) ({p^}i=i,2 are the densities of the subsystems) corresponding to a 
vibration of the superspace crystal in a direction between the average directions of the two 
sets of atomic surfaces. In the average, the atomic displacements of the subsystems are 
antiparallel and conserve the overall centre of mass. Of course, there are other possible 
modes corresponding to arbitrary angles 6. In composites we shall call sliding modes all 
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modes with r\ ^ (or equivalently 9 ^ 0, (3 ^ 0). Because r\ ^ 0, sliding modes involve 
a relative shift between the mass centres of the subsystems. The phason is a particular 
sliding mode which conserves the overall centre of mass, but there is a continuous set of 
other possible sliding modes for which the overall centre of mass is not conserved. 

The polarization parameter is zero for acoustical phonons (when rj = 0) and may in 
principle take any non-zero value for mixed phonon/phason modes. In the analytic regime, 
there are two hydrodynamic branches. The hermiticity of the dynamical matrix imposes the 
orthogonality of the modes belonging to the two branches, which reads: 



+A(3\K)(3 2 {K)<{^) 2 >=0 

for modulated crystals (33) 
{[1 +(3\K)][l + (3 2 (K)]+ < (^) 2 >} P W+ 

+{[1 - fi\K)][l - (3 2 (K)] + < (^) 2 >}p( 2 ) = 

for composites (34) 

In the non-analytic regime, the modulation functions are discontinuous and the phason 
branch has a gap (u p (K = 0) > 0). The acoustical phonon ((3 = 0) remains the only 
acoustical branch, and obviously no coupling is possible within the frequency gap (small K), 
because there is no phason there. Nonetheless, the coupling becomes possible for frequencies 
uj > uj p (0), leading to crossover phenomena (see Sec.f|). The derivatives of the modulation 
functions are singular. Hull functions may also become singular in this regime and the 



equations 22,23,24 should be understood in the sense of generalized functions. Although 



this issue deserves further study, it will not be addressed in this paper. 



3 Inner polarization and inelastic neutron scattering 

Lattice dynamics can be investigated by neutron inelastic scattering. 

In coherent neutron inelastic scattering, acoustical modes will be seen as branches issued 
from main and satellite Bragg reflections, whose positions in reciprocal space (for both 
modulated and composite crystals) are: 

fc r , s = 2ix\rja x + s/a 2 ] (35) 
Ignoring damping, the one-phonon inelastic coherent scattering cross-section obeys TO, M: 

(^fe) co ,~E P %-^))^yx 

x \Enbnexp[-W n (k)}( U P(K).k)ex P (-iky n )\ 2 (36) 

where K (already defined in Eqs. p0| , pl|) is the relative distance to a Bragg peak k TyS) 
i.e. k = k T ^ s + K. The index p is for different modes of angular frequency u p (K) and 
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displacements u^(K) in a normalized mode. b n is the average scattering length of the n th 
atom and exp[— W n (A;)] is the Debye- Waller factor. 

Let us consider that all scattering lengths and Debye- Waller factors are equal. Then, 
ignoring the frequency dependence l/u p (K) and the slowly varying factor k exp [— W n (A;)], 
the following quantity, that we call dynamical scattering factor, represents the contribution 
of one mode to the scattered amplitude along the dispersion branch to = u(K): 



1 N 

DSF(k) = lim —= > u n exp(—iky n ) (37) 
N ~*°° y/N i 

v n=l 



We have considered normal displacements Y2n=i \ u n\ 2 — !• The constant l/\/~N ensures 
the convergence of the sum. A pure acoustical phonon mode at k = has u n = I /y/N and 
DSF = 1. 

DSF contains information on both structure (via the atomic positions y n ) and dynamics 
(via the mode displacements u n ). 
The static structure factor is : 



1 N 

SF(k) = lim — } exp(-iky n ) (38) 

.W— >oo iv ^ 
n=l 

In Appendix 1 we have calculated the SF and the DSF close to a reflection fc ns , supposing 
that Eqs. |2^J23] , are accurate and that the modulation functions are smooth. It follows 
that for small modulation amplitudes (( = k r s sup |/| << 1) and for small reciprocal space 
distances to the reflection (% = \K/k r ^ s \ « 1), the ratio DSF/SF depends in lowest order 
of C, x only on the inner polarization, incommensurability ratio a and indices r,s (see also 

0)- 

For modulated crystals one has: 



\DSF(k r , s +K)\ 

\SF(k r ,s)\ 

1 

l+[l-t(A-)]2<(^.)2> 



ar+t(K)s 
ar+s 



-\ 2 



(39) 



where 



t(K) 



1-/3(K) 



1 + (3(K) 

The corresponding relationship for composite crystals is: 



(40) 



\DSF(k r . tS + K)\ 
\SF(k^)\ 



n\a + n 2 t 2 (K) 



ar + t(K)s 
ar + s 



-i 2 



(41) 
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where n%,n2 are the numbers of atoms of each subsystem within the periods ai,a,2, re- 
spectively (maximum values of indexes g, h in Eq. |3|). 

One should notice that in Eqs. j3i^ - |41| , |gpP depends on the interaction details only via 



the inner polarization parameter (3. For modulated crystals < (^-) 2 > also depends on the 
interaction details, but this quantity necessary for the normalization of the modes is the same 

I D S 

for all reflections k r ^ s . The K dependence of the ratio 1 1 is given by the K dependence of 
(3. This simple result can be used to interpret experimental data, but we should recall that 
our DSF differs from the experimental dynamical structure factor by a factor which depends 
slowly on k and that we ignored damping. 



4 Dynamics of composite crystals using the double chain 
model 



The double chain model was introduced [ITO, 15 1 in order to study the dynamics of composite 



structures made of two intermodulated subsystems f\. This model consists of two parallel 
chains of atoms. The atoms interact via pair potentials and move only longitudinally along 
the common direction of the chains. 

The Hamiltonian of the double chain is: 







„(2) (1) 
r m J &n ' 


y%) = 


tf } ) 2 

2m-! 


+ \ (y?> 


-J 1 ) _ 
Un-l 


ai) 2 ] + 


2m 2 


t(2) 

+ (y {2) 


(2) 
Vm—1 


a 2 ?] + 



(42) 



(1) (2) 

\^APm ) ' . , 2 \ (2) . s^-irfy™ ~~ y™ ■ 



where k^',k^ 2 ' and mi, m 2 are elastic constants and masses for the two chains, and V is 
the interchain potential of range r. 

4.1 Analytic regime 

The analytic regime corresponds to weak interaction between the subsystems, or equivalently 
to rigid subsystems (large k^ l \ k^ 2 '), and is characterized by continuous modulation functions 

The dynamical scattering factor for the double chain model with a Gaussian interchain 
potential 

(V(x) = —exp(—x 2 )) and equal masses mi = m 2 is shown in FigfJ. Several features may 
be noticed: 

• The dispersion curves of individual chains, touching the u = axis at the main Bragg 
reflections, represent the most significant features. 



2 A similar model was introduced in a different physical context M] 
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The intensity maxima form an hierarchical structure. 

From each reflection emerge two branches. One of the acoustical branches is strong at 
main Bragg reflections (r, 0) and weak at (0, s), while the other branch has the opposite 
behavior. At satellite positions ((r, s), r,s,^ 0), both branches are weak. 

For satellites, strongest DSF do not correspond to strongest SF. The strongest DSF is 
at (—1,4), while the strongest SF is at (1, 1). 

At high frequency there are nearly flat bands already observed for the Frenkel-Kontorova 



model or for the Fibonacci chain quasicrystal model [23 



In order to understand the low frequency aspect of Fig. |3] one has first to calculate the 
inner polarization parameter for the two acoustical branches using Eq.p7|. 

Fig.|| shows the dispersion curves and the inner polarization parameter (3 for the two 
acoustic-like branches. The acoustic-like branch 1 is the pure phonon (3 = 0, 6 = at 
K — and immediately it becomes a sliding mode concentrated on the first chain (/3 = 
1,6 — 72,7?.i = 1). At K — 0, the acoustic-like branch 2 is the phason with (3 = (p^ + 
j0 (i))/(p(2)_ p (i)) = _s±l ~ _ 4 .23, 6 = tan~ 1 [tan(-f 2 )(l+ P {1) / P {2) )] and at K ^ it becomes 
a sliding mode concentrated on the second chain because f3 = —1,6 = tt/2,1Z\ = 0. The 
displacements u n are shown in Fig.^. The jump of j3 very close to K = is a consequence 
of the fact that in our numerical simulations we use approximants and thus the phason gap 
is very small, but not strictly zero. Strictly speaking the inner polarization is undefined at 
K = because of the degeneracy, but in practice there is always a small gap that lifts this 
degeneracy. Furthermore, although the small K values where the acoustical phonon lives 
may be inaccessible to neutron inelastic scattering, the above remark is important when 
determining the speed of sound (this will be du/dK for K = 0,(3(K) = 0). 

The values of the inner polarization parameter explain why near main Bragg reflections, 
only one acoustic-like mode is strong (the one concentrated on the chain producing the Bragg 
reflection) and why the other is much weaker (see also Fig.^J). The mode with /3 — — 1 is 
visible near reflections r = and the mode with f3 — 1 is visible near reflections s = 0. This 
is expressed by Eq.^T] which has been illustrated in Fig.^. The above behavior of the inner 
polarization depends on the interchain potential and occurs when the slopes of individual 

dispersion curves of the two chains are not close one to another (v i = ^ = a\ yf^~ 7^ v 2 — 
= a 2\/^-)- The same inversion of mode intensities between (r, 0) and (0, s) main Bragg 



reflections can be found in Table. 1 of |24[ for the double chain model with a Lennard- Jones 
potential. 

When the single chain sound velocities V\, v 2 are close one to another 0, the inner polariza- 
tion behaves differently. Single chain modes couple strongly and tend to be non-concentrated, 
involving the participation of both chains (0 < TZ\ < 1 ) and thus contributing to DSF near 
both types of main Bragg reflections. This type of behavior is illustrated in Fig.|^ where we 
have chosen three different situations: v\ = 0.9t>2,fi = ^2,^1 = 1-1V2- Although the single 
chain sound velocities are slightly different in the three cases, after the interaction is set 



the meaning of "how close" depends on the type and range of the interchain potential. 
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in the dispersion curves look the same; this is the effect of mode repulsion. Nevertheless, 
the inner polarization parameters reach different values in the three situations. In the case 
v\ = t>2 one of the branches becomes very close to the pure phason (3 = —3.58, TZ\ = 0.34 
and the other branch is very close to the acoustical phonon (3 = 0.04, TZi = 0.65. In the 
other cases, there are two almost concentrated sliding modes: (3 = 2.23, TZ\ = 0.92 and 
(3 = -0.58,^i = 0.10 for v x = 0.9v 2 and (3 = 0.52, Tl x = 0.94 and (3 = -1.42, TL X = 0.05 for 
t>i = l.lv 2 . 



Once SF and (3 known, DSF can be calculated using Eq.[IT| The result is compared to the 
direct determination of the DSF (via Eq.[37]) for wave-vectors close to two satellite reflections 
in Fig.|5]. We have checked the validity of Eq.|41] close to all reflections with — 10 < r, s < 10 
and the result is illustrated in Fig.|[ The very good agreement is a proof that the undistorted 
atomic surfaces hypothesis is correct at least for large wavelengths. The errors in Eq.^T] have 
two origins: eventual deviations from the above hypothesis and the fact that the calculations 
leading to Eq.[Tl] are perturbative and valid for small modulation amplitudes. A direct test of 
the undistorted atomic surfaces hypothesis is to check Eqsf22|. [23|.|2|. This is illustrated in Fig. 
[T0| where we have plotted the ratio SSE/SST against K . SSE is the sum of squares error given 
by Eq.^, while SST is the total sum of squares estimating the variance of the hull function 
(SST = £ ft p/W(na a )- < >] 2 + EJU (2) (ma 2 )- < >] 2 ). R = yjl - SSE/SST 
represents the multiple correlation coefficient (R close to 1 means good quality of the linear 
regression [pOfl). 



4.2 Non-analytic regime 

The non-analytic regime occurs when the interaction between the subsystems is strong, or 
equivalently when the subsystems are soft (small k^), and is characterized by discon- 
tinuous modulation functions jn|. As a consequence, the phason branch starts with non-zero 
minimum frequency (gap at K — 0). This gap can not be noticed in Fig.|3| (it is too small), 
but it can be seen in Fig.^. Several other gaps exist at higher frequencies and K ^ 0. 

The variation of (3 inside the two branches is shown in Fig.|]. For low frequencies the 
branch 1 has (3 = 0, hence it is the pure acoustical phonon, while the branch 2 is the phason 
[(3 = (p( 2 ) + p(!))/(p( 2 ) - p (1) )). This can also be seen in Fig.[7| at K = 0. The branch 
1 corresponds to uniform displacements of the atoms, the branch 2 has the rapid phason 
oscillations and a sliding character. For higher frequencies both branches are sliding modes, 
almost concentrated on single chains. Like in the analytic regime, mode concentration can 
be avoided when the slopes of the individual dispersion curves are similar. 

The difference with respect to the analytic case is that the interval of values of K over 
which the inner polarization changes is much larger (it scales like the gap). A change of 
the slope du/dK accompanies the change of (3 as can be seen in FigJJ. This could have 
consequences when different methods of investigation are employed. The speed of sound 
determined by ultrasound and light Brillouin scattering in the K = limit measurements 
may be different from the slopes du/dK determined by neutron measurements. 

Although the static structure factor is different from that in the analytic regime, the 
ratio \DSF/ SF\ is still given with reasonable accuracy by Eq.41. This is shown in Fig|| 



where the numerical DSF, and the DSF calculated from SF via the Eq.[y] are compared 
for two satellite reflections. The accuracy of Eq.|4l] was tested for all reflections (r, s) with 
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— 10 < r, s < 10. The result is presented in Fig.[]. 

To summarize, in both analytic and non-analytic regimes, low frequency features (inten- 
sity of acoustical branches close to different reflections) of DSF depend on the SF and on 
the value of the polarization parameter (3. Both SF and Q depend on the interaction details 
of the system, furthermore may also depend on frequency. For various model parameters 
and K values, (3 scans a broad domain of values, corresponding to inner polarizations of 
various inclinations 9 within (— tt/2, 7r/2]. 



5 Dynamics of modulated crystals via the DIFFOUR 
model 

DIFFOUR(discrete frustrated $ 4 ) models simulate |25| displacively modulated crystals with 
competing interactions. A simple example is represented by a chain of atoms in positions y n 
such that the potential energy depends non-linearly on the fluctuations x n = y n — y n -i — a\ 
of the first neighbor distances. The Hamiltonian is: 



H(p n , Vn) = 

Pi , A „2 , 1 4 n„ „ , n„„ \ ( 4 3) 



E 



For A > A* = B 2 /(AC) + 2C the ground state is not modulated (x n = 0). Modulation 
occurs for A < A* with a periodicity a 2 that can be incommensurate to the periodicity a\ 
of the basic structure: a 2 = 2nai/Q, where Q = arccos(^j). Close to the commensurate- 
incommensurate (C-IC) transition the modulation is analytic (A w A*). Analyticity breaks 
(the modulation function becomes discontinuous) when A becomes too small. For A ^ A*, 
the slopes of the hydrodynamic branches are given by = { — [(2B + A + C) cos K — 4(2C + 
-B) cos2.fr + 9C cos 3K]} -1 / 2 , with K = for the acoustical phonon and K = Q for the 

phason. These slopes are equal (4tp) , = (4^) , = 3<i/— when B = —2C. 

1 1 1 \ aK ) phonon \ dK ) phason y m 



In Fig.rTT] we have represented the dispersion curves and the inner polarization parameter 
8(K) in three situations: , ~ i^w) , for an analytic modulation, (4^) , ^ 

1 y ' V dK I phonon \ dK I phason u 1 \ dK I phonon ' 

(-757) l for an analytic and a non-analytic modulation. In all these situations 3 is constant 

\ aK ) phason J J ' 

and equal to for the acoustical phonon and is equal to —1 for the phason. No visible 



phonon/phason coupling occurs. Further insight is given by the hull functions (Fig. O) 



With very good accuracy the phason hull function is the derivative of the modulation function 
and has thus Q — — 1. The acoustical phonon corresponds to uniform displacements of the 
atoms. A small amplitude modulation (though with a wavelength that is a fraction of a 2 , 
hence different from the derivative of the modulation function) add to them when phonon 
and phason dispersion curves are close (Fig. |T2"| a)). 

The DIFFOUR model, like the DCM model fit to our general scheme and satisfies the 
non-distorted atomic surfaces condition. Contrary to the DCM model, the values of the 
inner polarization parameters are practically restricted to Q\ — and j3 2 = — 1, correspond- 
ing to the acoustical phonon and to the phason. The following simple argument suggests 
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that close to the C-IC transition the phason/phonon coupling vanishes not only for the DIF- 
FOUR model but for modulated crystals in general. Considering that the potential energy 
1S Emn^fen — that c l° se to the transition the modulation is sinusoidal and of small 
amplitude, that phonon displacements are u n ~ exp(iknai) and that phason displacements 
are v n ~ f (nax)exp(iknai) ~ exp[i(k + 2ir / a2)nai], then the phason/phonon coupling scales 

like Y,m,n U n V "( na l ~ ~ Y,m exp(2inma 1 / d 2 ) = 0. 



6 Composites vs. modulated crystals 

Is it possible to distinguish between composites and modulated crystals by studying their 
dynamics? A possible experimental procedure would be to scan DSF in the neighborhood 
of the main or satellite reflections. One may choose a set of reflections and determine the 
ratio \DSF/SF\ for the acoustical modes. From Eqs. [M O it follows that the dependence 



of the \DSF/ SF\ ratio on the reflection indexes r, s has the same form for modulated and 
composite crystals. The only way to distinguish between different classes of aperiodic crystals 
is via the polarization parameters. We have shown that these are different for modulated 
and composite crystals. 



First of all, the orthogonality relations 33,34 that relate inner polarizations of the acous- 



tical modes read differently for the two classes. For instance a j3i — — 1,6*1 = tt/2 mode 
(phason in modulated crystals and concentrated sliding mode in composites) is accompanied 
by a f3 2 = 0, 6*2 = (acoustical phonon) in modulated crystals and by a /3 2 — 1)#2 = 72 
mode (concentrated sliding mode) in composites. The (3 2 = and (3 2 = 1 modes scatter very 
differently and this should be seen in neutron inelastic scattering experiments. 

On the other hand, close to the C-IC transition, the inner polarization of normal modes 
in modulated crystals is restricted to /3i = —1 (phason) and [3 2 = (acoustical phonon). 
Although the numerical simulations using the DIFFOUR model showed only negligible pha- 
son/phonon coupling, the theoretical arguments that we employed do not exclude the pos- 
sibility of a phason/phonon coupling far from the C-IC transition, for instance in the non- 
analytic regime. In composites, various inner polarizations are possible. 

Let us see how the above theoretical considerations can be used to interpret available 
experimental results on structure and dynamics of incommensurate modulated and composite 
crystals. 

Incommensurate modulated crystals have been studied rather intensively and there are 
some well characterized model systems |26, 27, 28, 2£]. 



For the incommensurate compound ThBr^ |3(j] the modulation function is analytic, al- 
most sinusoidal (only first order satellites are visible). In this case there is no difficulty 
to identify the main reflections and the classification task is trivial. A gapless phason was 
observed by inelastic neutron scattering close to the first order satellite (2, 3, 0.69) of the 
main Bragg (2,3, 1). Despite that the intensity of the studied satellite reflection is strong, 
the acoustical phonon branches were measured only around the main reflections. 

In the phase III of biphenyl |2(| one can observe strong satellites, up to order 3 in 



neutron scattering, suggesting a non-analytic (soliton) regime or (which is equivalent) a 
higher-order commensurate phase. Coherent inelastic neutron scattering scans around a 
first-order satellite have revealed a frequency gap in the phason branch at K 7^ 0. The 
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presence of a K = gap that would be expected in the non-analytic regime has not been 
detected, perhaps because of the instrumental resolution and the heavy damping of these 
excitations. 

There are many structures that can be named " composites" . Still, a precise classification 
is not established yet. 

The best studied composite system is the mercury chains compound Hg 3 _sAsF 6 with 
8 ~ 0.18 [32|]. In this case there is no doubt that the structure is composite, being made 
of two weakly interacting subsystems. The tetrahedra of AsF e form the host tetragonal 
subsystem. There are two mercury subsystems in two orthogonal sets of channels of the 
host, running along the equivalent incommensurate a and b axes of the tetragonal system. At 
room temperature the host main Bragg peaks are easily recognizable at positions (h, k, I) and 
mercury chains produce diffuse thin sheets at (n(3 — 5), k, I) orthogonal to a* and at (h, n(3 — 
S), I) orthogonal to b*. Near main Bragg peaks of the host and along the incommensurability 
direction a, acoustical longitudinal modes with a speed 4& = 2130ms" 1 were measured. Near 
mercury chain diffuse sheets and along the same direction a gapless longitudinal branch was 
measured with a very different speed ^ = 3616ms" 1 . According to our discussion this 
suggests the existence of two acoustic-like modes with inner polarization parameters fti — 1 
and fa — ~ 1) concentrated on the host and on the mercury chains, respectively. This picture 
is consistent with the weak inter-subsystem coupling and the large difference between slopes 
of the dispersion curves (Figs.f|,[5|). 

The alkane/urea inclusion compounds |33| are also guest/host type composites with two 
easily recognizable subsystems. At all studied temperatures, guest chains (alkane) have only 
ID order and produce diffuse sheets, like in the case of mercury chains, but the alkane 
modulation is found to be strongly non-harmonic |36|, |37| . There are no complete studies 
of dynamics for this compound, nevertheless special features of the longitudinal phonons 
along the incommensurability direction were observed near main Bragg reflections of the 
host by neutron scattering |34| and confirmed by Brillouin scattering f38fl . These phonons 
have unusually large damping (approx. 400 GHz, independent of K ||35|| ) compared to the 
other acoustical phonons (damping scaling like K 2 ) that could be explained by a sliding 
character (finite /3, phason/phonon coupling). Although this was not studied in this paper, 
some authors |L3], [L4| suggest that damping and sliding are connected via internal friction. 
This internal friction can have an intrinsic origin (occurring for large relative velocities |39] ) 
or an extrinsic one (due to defects). 

The case of Bi — 2212 lamellar superconductor is more complex because there is no 
clear separation of the subsystems and the structural classification is still controversial 
[0, f40|. Neutron inelastic scattering f| showed the presence of two longitudinal acousti- 
cal branches close to two strong reflections and along the incommensurability direction. 
The slopes of the dispersion curves are ^ = 2400 and 5900ms -1 and their intensities do 
not change when one passes from one reflection to another. According to Eqs. [39|J40|r4I 



(DSF/ SF) 2 Q J \DSF/ SF) 2 rQ ^ = t 2 irrespective of the type of structure, hence t 2 should be 
close to 1 for both acoustical branches. This is possible only if << 1 and l/^l >> 1- One 
could not speak of a modulated structure, because this contradicts Eq.|33|. If one of the modes 
is the acoustical phonon {(3i — 0, t\ = 0) then in modulated crystals the other mode must 
be the phason (/? 2 = — 1,^2 = oo) being thus invisible close to main Bragg peaks (s = 0). 
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In Ref. [0] it was assumed that the structure is composite and that the two reflections are 
the main Bragg peaks (2,0) and (0, 1). The intensity ratios are compatible with Eq.[34] and 
can be explained by considering that one of the longitudinal modes is the acoustical phonon 
(Pi = 0) and that the second is the phason (\p 2 \ = usjyi = 134.33 >> 1 (this was computed 
considering that the densities for the subsystems Sr<iCaCuiO% and Bi 2 02.2\ are 3.248g/cm 3 
and 3.2g/cm 3 respectively). For both the above inner polarizations the DSF/SF ratio is not 
changing between the r = and the s = Bragg reflections (Fig.|B|b)). To conclude, dynamics 
suggests that the superconductor is composite, but the modes are non-concentrated, which 
is different from the case of mercury chains compounds. This situation is reproduced by the 
DCM model in the rather special situation when v\ ~ v% (Fig.|8|). Before deciding which is 
really the case, further study is needed in order to check how the interaction details (range 
and type of the inter-chain potential) and the dimensionality (DCM generalizations to 2D, at 
least) influence the inner polarization. Referring now to the regime, some of the modulation 
functions used in the structure refinement have large amplitudes and are non-analytic (in 
particular an occupational modulation is introduced for the oxygen atoms in the BiO layer 
of the second subsystem) , while other modulation functions could be analytic. One may thus 
think of a situation when the phason gap is below the experimental resolution of 10 GHz. 
As discussed in Sec. f|, this hypothesis could also explain the observed discrepancy between 
the ultrasound measurement of the speed of sound (4200ms" 1 ) and the slope of the neutron 
dispersion curves. 

The dynamics of spin ladder compounds (Sri4_ x Ca x )Cu24,0 4: i [fll|] has been very recently 
investigated ||42|| . Again two low frequency longitudinal modes (along the incommensurabil- 
ity direction)were observed close to a main Bragg peak of each subsystem Cu0 2 (chains) 
and (Sr,Ca)2Cu20 3 (ladders). One mode is acoustic-like (no gap), the other has a gap of 
0.35THz (for x — 0) and can propagate also transversally. Interestingly enough, when the 
latter propagates longitudinally it looses intensity very rapidly with K. Because of the gap, 
the mode at low frequencies should be the acoustical phonon. The loss of intensity with 
K of the second mode (the phason) suggests that this one tends to become a sliding mode 
concentrated on one of the subsystem when K is increasing, like in Fig.||]. The orthogo- 
nality conditions impose that the acoustical phonon should also change character when K 
increases, concentrating on the other subsystem. Our theory predicts that the observation 
(not yet performed) of a Bragg peak belonging to the other subsystem could reveal an op- 
posite behavior of intensity of the acoustical phonon (this will increase with K, see Fig.|5|, 
non-analytic regime, (1,0) and (0,1) main Bragg peaks). Of course, a more cautious analysis 
should also include damping effects. 

7 Conclusion and Discussion 

We emphasized the utility of the concept of inner polarization in the description of low 
frequency excitations in composite and in modulated aperiodic crystals. This concept has 
a meaning for long wavelengths in the analytic and even weakly non-analytic dynamical 
regime, generally for small modulation amplitudes. The inner polarization is related to 
internal degrees of freedom of the aperiodic crystals and should not be looked at as the 
direction of atom displacements (this is the physical polarization). Modes having different 
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inner polarizations differ by the amplitude of rapid oscillations added to a uniform overall 
displacement in modulated crystals and also by the "sliding" character and the relative 
participation of the subsystems in composites. The inner polarization can be defined for 
modes that do not change the shape of the superspace atomic surfaces. The possible values 
of the inner polarization parameter (3 are different in composites and modulated crystals. 

In modulated crystals, close to the C-IC transition the phason and the acoustical phonon 
should be uncoupled. 

In composites, the phason and the acoustical phonon can be coupled. We called "sliding- 
modes" all the modes with (3 ^ (the only exception is the acoustical phonon) because they 
produce a relative shift of the mass centres of the subsystems. The pure phason is a special 
sliding mode, the only one that conserves the overall mass centre. One can distinguish two 
extreme types of dynamics in composites: concentrated modes (two sliding modes with (3 G 
{ — 1,1}, involving each one displacements of a single subsystem) and non-concentrated modes 
(the acoustical phonon, the pure phason or other sliding modes involving both subsystems). 
Crossover is possible from one type of dynamics to another. The variation of the inner 
polarization and of the participation ratio with the wave vector K is a different phenomenon 
from the crossover between diffusive and propagating modes discussed in [j43[ , |T3[ fL4|j . 

The inner polarization of the modes influences the neutron inelastic scattering. For 
modulated structures the pure phonon is visible near all main Bragg reflection and satellites, 
but the phason is visible only near satellites. For composites, the "concentrated" modes are 
visible only near main Bragg reflections of the subsystem on which they are concentrated, 
or near satellites, "non-concentrated" modes are visible near both types of main Bragg 
reflections and also near satellites. 

The above properties can be used to distinguish between modulated structures and com- 
posites and also between the two extreme types of dynamics in composites. More information 
could be (at least in principle) extracted by an investigation of the satellites. The inner po- 
larization can be determined from the ratio between the dynamical structure factor and the 
statical structure factor at various reflections. Of course, for more precision one has to con- 
sider also the neglected factor kexp[—W(k)\ in the DSF and the effect of damping on the 
line shape. 

In experiments,it is important to know for a given mode the set of satellites where this 
is visible. A pure phonon scatters strongly near strong satellites. For modulated structures 
the empirical observation is that the same thing is true for the phason. This can by justified 
by multiplying eq.|39| for (3 = —l,t = oo and k 2 exp[—2W (k)} where k ~ (rat + s). We 
get I ~ s 2 \SF\ 2 exp[— 2W(k)}, where / is the intensity scattered by the phason. If \SF\ 2 
decreases with s quicker than 1/s 2 , then the first order satellites s = 1 are the most intense 
and correspond to the strongest phason inelastic scattering. For composites, from eq.^1] we 
get: 

• for (3 = —1 (sliding mode concentrated on the second subsystem) I ~ s 2 \SF\ 2 exp[— 2W{k)]. 

• for (3=1 (sliding mode concentrated on the first subsystem) I ~ r 2 \SF\ 2 exp[— 2W(k)}. 

• for (3 = p( 2 )^p(i) (pure phason, hence non-concentrated) / ~ (r— -^^s) 2 \SF\ 2 exp[— 2W(k)}. 

We face here a more complex situation. In composites, one can have significant values of 
the SF for high values of r or s for satellites of the type (1, s) or (r, 1) (which are first order 
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for the second and the first subsystem, respectively). Even if the (1,1) satellite could be 
statically more intense, the scattered intensity may be stronger for (1, s) (for f3 = —1) or for 
(r, 1) (for (3=1). Furthermore, for non-concentrated sliding modes, there is an extinction 
rule: the mode will be invisible close to satellites satisfying r = -^-j^-s (in the DCM model 
with mi = rri2 this becomes r = s). This type of extinction has been previously reported in 



theoretical models [ I9| . Thus, sliding modes may not give the strongest inelastic scattering 
near the strongest satellites. In a given practical situation one should use Eq.|41| together 
with the measured values of \SF\ in order to calculate \DSF\ and to find out where this is 
the strongest. 

From our discussion of the inner polarization it becomes clear that sentences of the type 
"this is the phonon mode and this is the phason (or sliding) mode" referring to acoustic- 
like modes in composites are at least incomplete. Sliding modes may have various inner 
polarizations and the pure phonon may simply not exist when there are two sliding modes. 

Further study is needed in order to extend (or to specify the limits of) the validity of the 
results presented here, including the cases of large amplitude and occupational modulations. 
For a better analysis of the experimental results a treatment of damping is necessary. 

8 Appendix 

Let us define the following functions : 

^k^) = ^M-ikf^(x)} (44) 

d^*) = ^0*) ex PH<3(^)] (45) 
rfkg and 4>^ g are periodic of period a 2 , while ipjp h and <jf£\ are periodic of period a\. Let 

k,(g,h),m an< ^ ^k,(g,h),m 



^kfJh) m an d $kfah\ m ^ e the Fourier coefficients of the above functions, such as: 



^k$,h)( x ) = E m ^(^ ) , m exp(27rzmx/a {2> i ) ) (46) 
^ { k,S,h)( x ) = E m ^(g!/ l) , m exp(27rima;/a( 2 ,i ) ) (47) 

Let us consider periodic approximants such that a 2 /ai = p/q, p,q relatively prime inte- 
gers. The final result will not depend on this periodic boundary condition because at the 
end we shall impose p, q — > oo, a 2 /ai — > a. 

In composites, the zero frequency normalized acoustical modes are: 

*£U = N m[! + + ^^fK + 4 1} )]x (48) 
x exp[iK(nai + Sg 1 ^)] 

u ( Z,k = N p^ ~ /W ~ W{K) d ^{ma 2 + 4 2) )]x 

xexp[iir(ma 2 + 4 2) )] (49) 
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where 



= wtt 1 + ^ K ')f + 4 ^ 2 (#) < l^l 2 >}+ 

+n 2 g{[l - /3(K)] 2 + 4/3 2 (iT) < |^| 2 >]} 



ni,n2 being the number of atoms per each sub-period. From Eq.3? and k = 2-n{r/ai + s/a2) + 
K one gets: 



a i 



9=1 



+ 2/w*S,j + EK 1 - - 2W*S, r ]} 



(50) 



where 



p,q~^oo q 

,Ai) 

(ma + n 2 ){rna[(l + /?) 2 + 4/? 2 < (^-) 2 >} + (51) 



+ n 2 [(l-/3) 2 + 4/3 2 <(^H 2 >]} 

Similarly, we obtain: 



I "1 "2 

^ raffiP *(M = ^[«EC+E€,J (52) 

9=1 h=l 

In order to obtain the above results we have used the following version of Poisson's sum 
rule: 



1 P 

lim - > f(naq/p)exp(—2irinsq/p) = f s (53) 

p^oo r> — ' 
K n=l 

that is valid for any continuous function /, periodic of period a. In order to prove Eq.[53| 
one could use f(x) = f s exp(2iT i sx/a) and Yln=i exp[27rinmq/p\ = pJ2k^"i,k P , valid for 
m, p, q integers p, q relatively prime. 

For small modulation amplitudes |/c/ ( - 1 ' 2 ' ) | << 1, we can use the following linear approx- 



imations of Eqs.|44,|45 
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fL « 5r,o - S r » 2^/^. Substituting in Eqs.|0f2] one gets : 



DSF composites (k) = N- l {a,f3{K)){an x {l + ( 5(i^)]5 S)0 + 
+n 2 [l - /Wl^o - iK/3(K)F- s>a - 
-»[27r(l + + 2tt(1 - + if]F+ Si J (54) 

^composites = ..^[^^ + 

where 



ni ri2 

= « E © ± E S ( 56 ) 

9=1 /i=l 

The same kind of arguments were used for modulated crystals to obtain: 



Demodulated^) = [\ e (K)\ 2 + \r]{K)\ 2 < (J) 2 >]-^ 2 X 

x{e(K)5 sfi -t[e(K)(K + 27r^) + 
(e(K) - V (K))2n±] < f 9jS >} (57) 

^modulatcd^) = _ 2m( j_ + ^ < f gg > (5g) 

where < f 9jS >= ^ YT g U fg,s- 
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Figure 1: Superspace embedding of a) modulated incommensurate crystals (one set of 
atomic surfaces) b) composite incommensurate crystals (two sets of atomic surfaces). For 
illustration, we have chosen sinusoidal modulation functions. 6 is the angle made by the 
inner polarization and the physical line (horizontal). The inner polarization parameter is 
(3 = tan(Q) /fttan^) — tan{6)}. The acoustical phonon corresponds to 6 = 0,f3 = 0. For 
composites, modes concentrated on the first subsystem have inner polarization parallel to 
the atomic surfaces of the second subsystem, i.e. 9 = 72,/? = 1. Modes concentrated on the 
second subsystem have inner polarization along the atomic surfaces of the first subsystem, 
i.e. 6 = tt/2,/3= -1. 




Figure 2: Different inner polarizations P and displacement hull function U for modulated 
crystals; a) phason with P _L U, b) acoustical phonon P||U (U is constant), c) mixed 
phonon/phason, d) distorted atomic surface, inner polarization can not be defined. 
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Figure 3: a) DSF for the double chain model (a = (l + v / 5)/2,mi = m2,^ ! V^ = 1.2, short 
range Gaussian potential r = ai/3). The calculations were performed on approximants of 
period a = pai = qci2, p = 233, q = 144. Pseudo-colors show DSF values in a logarithmic 
scale, b) DSF for the two acoustical modes are compared to SF at different reflections, uo is 
in y/k^> /mi units 
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k^'/k^ 2 ' = 1.2, short range gaussian potential r = ai/3). a) Folded dispersion curves within 
the first Brillouin zone of the approximant of period a = pa\ = qa2, p = 233, q = 144. u 
is in a//^ 1 ) /mi units; b) Inner polarization parameter j3 along the two branches. Analytic 
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Figure 5: Numerical values of DSF near the satellites and main Bragg reflections, along 
the low frequency branches, as a function of the distance in the reciprocal space relative to 
the reflection. Thick lines were obtained using SF and Eqj41|. Same parameters as for Fig.|j. 
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Figure 6: DSF/SF ratio for the main Bragg reflections of composites as a function of the 
inner polarization parameter of the mode, a) for a = (1 + v / 5)/2, as in Fig.[|; b) for a = 2.2 
as in the Bi — 2212 lamellar superconductor. 




Figure 7: Lowest frequency modes in the double chain model. Same parameters as for Fig.|j. 
The displacements are longitudinal and correspond to n = 1, 233 for the first chain and 
n = 234, 377 for the second chain. They are uniform for the acoustical phonon, antiparallel 
for the phason and practically involve only one of the two chains for concentrated sliding 
modes (the other chain presents only the small amplitude rapid oscillations equal to the 
derivative of the modulation function). 
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Figure 8: Unfolded dispersion curves and polarization parameters in the DCM model, 
when the single chain sound velocities Vi = Oiy/k^/rrii, % = 1, 2 are close one to another, a) 
vi = 0.9v 2 ; b) v\ = v 2 ; c) t>i = l.lv 2 ; 
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Figure 9: DCM model: DSF obtained directly from Eq.|37| vs. DSF calculated from SF via 
Eq.[|l] for various main and satellite reflections. 
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Figure 10: DCM model: Direct test of Eqs. |23|j24] . The ratio SSE/SST for the two acoustical 
branches in the analytic regime is plotted against K. The quality of the regression is excellent 
(SSE/SST « 1 ). 
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Figure 11: DIFFOUR model: Unfolded dispersion curves and polarization parameters. We 
used periodic boundary conditions pai = qa 2 . a) C = 0.4, A = 0.99A*, B = —2.1C, 
p = 50, q = 17, analytic regime, b) C = 0.4, A = 0.99A*, B = 2.5C, p = 50, q = 7, analytic 
regime, c) C = 0.4, A = 0.65 A*, B = 2.5C, p = 50, q = 7, non-analytic regime. 
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Figure 12: DIFFOUR model: real and imaginary parts of the actual (*) and predicted (x) 
(by Eq.|22l) hull functions for the phonon (up) and for the phason (down) at K = 0.8 — . a) 
C = 0.4, A = 0.99A*, B = -2.1(7, p = 50, q = 17 analytic regime, b) C = 0.4, A = 0%9A*, 
B = 2.5(7, p = 50, q = 7 analytic regime. 
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